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Abstract 

We investigate the liquid-gas second-order phase transition in liquid ^He using the variational 
calculations based on the cluster expansion of the energy functional. We also compute the critical 
point exponents of liquid ^He which are in agreement with experimental data. 



I. INTRODUCTION 



The liquids ^He and ^ife are the only quantum liquids which exist naturally. The word 
"quantum liquid" comes from the fact that for these systems, the interatomic distance is at 
the order of their De Broglie wavelenght. ^He is a liquid of strongly interacting fermionic 
atoms which behaves quite differently than the normal liquids at low temperature (Wilks, 
1970; Kent, 1993). The properties of liquid ^He have been studied using different many- 
body techniques (Clark and Westhaus, 1966; Nafari and Doroudi, 1995; Luijen and Meyer, 
2000; Kindermann and Wetterich, 2001; Takano and Yamada, 1994; Viviani et al., 1988; 
Pricaupenko and Treiner, 1995; Fantoni et al., 1982; Krotscheck and Smith, 1983; Priman 
and Krotscheck, 1982). Recently, the behavior of hquid ^He near its critical point has been 
investigated using path-integral molecular dynamics and quantum virial expansion (Miier 
and Luijten, 2002). 

One of the most powerful techniques in many-body calculations is the variational method 
which is based on the cluster expansion of the energy functional (lowest order constrained 
variational method) (Owen et al., 1977; Bordbar and Modarres, 1997, 1998; Modarres and 
Bordbar, 1998; Bordbar and Riazi, 2001, 2002; Bordbar, 2002a, 2002b, 2003, 2004; Bord- 
bar and Hashemi, 2002). This is a fully self-consistent method and does not introduce any 
free parameter to the calculations. The crucial point in this method is the functional min- 
imization with respect to the two-body correlation function subjected to the normalization 
constraint which finally leads to a Euler-Lagrange differential equation. The convergence of 
its results has been shown by computing the three-body cluster energy term (Bordbar and 
Modarres, 1997). 

The liquid-gas phase transition near the critical point (second-order phase transition) is 
an interesting subject in statistical mechanics. This behavior, critical phenomena, is caused 
by the existence of singularity in thermodynamic functions of the system at the transition 
point. The nature of these singularities in various measurable quantities at the critical point 
is described by the critical exponents. 

In our previous paper, we have calculated some thermodynamic properties of liquid ^He 
using the variational method shows a nice agreement between experimental data and calcu- 
lated results especially for free energy and entropy (Bordbar and Hashemi, 2002). In this 
article, we present the critical behavior of the liquid ^He. We organize the paper as fol- 
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TABLE I: Critical point properties of liquid ^He. 





Tc{K) 


Pc{A-^) 




Our results 


4.36 


0.0054 


0.0139 


Exp. results (Heller, 1967) 


3.324 


0.00834 


0.00844 


Exp. results (Pittman et al., 1979) 


3.317 


0.00827 


0.00846 



lows: In section [TTl we obtain the critical properties of liquid ^He by calculating the critical 
isotherm. We investigate the critical behavior of liquid ^He by computing the critical point 
exponents in section IIIII 



II. CRITICAL ISOTHERMAL EQUATION OF STATE 

The equation of state is the key point for investigating the second-order phase transition in 
a hydrostatic system. The isothermal equation of state can be calculated from the Helmholtz 
free energy, F: 

p^pfl., (1) 

where P, T and p = y are the pressure, temperature and number density, respectively. 
and V are the total number of particles and volume. To obtain the free energy of liquid 
^He, we use the variational method explained in Appendix. The result for the equation of 
state at the critical temperature (critical isotherm) is shown in Fig. ([T]). As seen from the 
Fig. ([T]), at the critical point, the isotherm curve shows an inflection point which satisfies: 

dP. d'^P. 



dp (9p2 

where Tc is the critical temperature. The calculated critical temperature, density (pc) and 
pressure (Pc) of liquid ^He are presented in Table [B The experimental results (Heller, 
1967; Fisher, 1967; Pittman et al., 1979) are also given for comparison. We can see a good 
agreement between these results. 
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FIG. 1: The critical equation of state for liquid ^He. 
III. CRITICAL EXPONENTS 

For a hydrostatic system, the two-phase coexistence conditions are 

p p 

liquid ^ gas 
IJ'liquid fJ'gas ) 

where the Huquid and ^gas are the chemical potential of liquid and gas phases respectively. 
As the temperature increases, the liquid density decreases and the gas density increases. At 
the critical temperature these densities become equal to each other. This behavior for ^He 
is shown in Fig. The order parameter puquid — Pgas which is defined to investigate the 
critical behavior of this system vanishes at the critical point. However other thermodynamic 
properties diverge at this point. The critical point exponents are defined to study the 
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FIG. 2: The liquid and gas densities versus temperature for ^He. 

asymptotic behavior of singular thermodynamic functions near the critical point. For this 
purpose, the following functions for the thermodynamic quantities are introduced (Garrod, 
1995): 

• Order parameter 

We can define the exponent /? for this parameter as follows: 

Pliquid - Pgas ~ {-^Y] ^ ^ 0~, (4) 

where 

(5) 

The critical exponent {3 characterizes the behavior of the order parameter and of 
course, the above function is meaningful only below the critical point in the region 
where the order parameter is not zero. To obtain /?, we draw the order parameter as 
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FIG. 3: The order parameter versus e on log-log scale for ^He. 

a function of e on the log-log scale in Fig. ([3]). The slope of this figure yields the value 
of /3 = 0.56239 ±0.01386. 

• Pressure 

By defining the exponent (5, we can describe the critical isotherm 

F-P,~(p-Pe)^; p-^pc, (6) 

where e = (T = T^. In Fig. (jlj) , F — Pc as a function oi p — pc is shown. The value 
of 5 obtained from this figure is 3.31032 ± 0.08192. 

• Heat Capacity 

The exponent a' and a characterize the behavior of specific heat (Cy) below and above 
the critical temperature respectively along the critical isochore (V = Vc) 
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FIG. 4: The P — Pc versus p — pc sX critical temperature (Tc) for ^He. 

Cy. = (6)-"; 6-^0+ (7) 

In Fig. ([5]), the specific heat along the critical isochore versus e is shown. The values 
a' = 0.1018 ± 0.0001 and a = 0.10609 ± 0.0014 are extracted from the Fig. ([5]). 

• Isothermal Compressibility 

For describing the behavior of isothermal compressibility {K) near the critical point, 
the exponents 7 and 7' are defined to be 

K = (e)-^; e 0+. (8) 

The calculated values of isothermal compressibility shown in Fig. ([6]) leads to 7 = 
1.05343 ± 0.01077 and 7' = 1.05343 ± 0.01077. 

We have presented the whole critical exponents for the ^He in Table HTl The experimental 
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FIG. 5: Specific heat along the critical isochore as a function of e above (full curve) and below 
(dashed curve) critical temperature for ^He. 



TABLE II: Critical exponents for ^ife. 





& 


b 


a' 


a 
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Our results 


0.5624 


3.3103 


0.1018 


0.1061 


1.0534 


1.0560 




±0.0139 


±0.0819 


±0.0001 


±0.0014 


±0.0108 


±0.0093 


Exp. results (Heller, 1967) 


- 0.361 


4.21 


~ 0.105 


~ 0.105 


~ 1.17 


~ 1.17 


Exp. results (Pittman et al., 1979) 


0.322 ± 0.002 








1.19 ±0.01 





results (Heller, 1967; Fisher, 1967; Pittman et al., 1979) are also given for the comparison 
in Table [Tll There is a good agreement between our calculations for the critical exponents 
and the experimental results. From Table [Tll it can be seen that the Griffiths and Rush- 
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FIG. 6: Isothermal compressibility as a function of e above (full curve) and below (dashed curve) 
critical temperature for ^Re. Ki is the ideal fermi gas compressibility at p = Pc and T = Tc. 

brooke inequalities (Huang, 1987; Griffiths, 1965) are satisfied by our results for the critical 
exponents of ^i?e, 

a + 2/3 + 7 > 2 
a + /3(l + (5) > 2. (9) 

IV. SUMMARY AND CONCLUSION 

The liquid-gas phase transition near the critical point is of special interest in statistical 
mechanics. In this work, we have computed the critical equation of state for liquid ^He which 
led to critical density, temperature and pressure of this system. The critical exponents, 5, 
a and 7 for this system are computed. The calculated critical exponents satisfies the Griffiths 
and Rushbrooke inequalities. A comparison between our results and experimental data 
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is made which shows a good agreement between theoretical calculation and experimental 
results. 
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APPENDIX 

In this appendix, we give a brief review to obtain the free energy of liquid ^He using the 
lowest order constrained variational method based on the cluster expansion of the energy 
functional (Owen et al., 1977; Bordbar and Modarres, 1997, 1998; Modarres and Bordbar, 
1998; Bordbar and Riazi, 2001, 2002; Bordbar, 2002a, 2002b, 2003, 2004; Bordbar and 
Hashemi, 2002). In this method. We choose a trial many-body wavefunction as 

i<j 

where f{ij) is the two-body correlation function and $ is the Slater determinant of nonin- 
teracting particles wave-functions (plane waves). We then apply the cluster expansion to 
the energy per particle (Clark, 1979) and keep one and two-body energy terms, 

£; = ' ^ (11) 

where 

Ei = E^n{k,), (12) 

i 

E2 = ^E<^^'I«'(12)|U-J^>- (13) 
In the above equations, n{ki) is the Fermi-Dirac distribution function and 

w{l2) = ^(Vuf{l2)y + f{12)V{12), (14) 

where T^(12) is the interatomic potential. In the thermodynamic limit, Eqs. (IT^ and (I13p 
read: 

i-oo 

^1 = [ nik)kUk, (15) 
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where p is the number density, 
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+ -,f{r)V{r) 



1 r°° 
p = — n{k)k^dk, 
vr^ Jo 



1 /7(r) 



r dr, 



7(0 = A / 

TT^ JO 



°° sin kr 



kr 



n(k)k^dk, 



(16) 

(17) 
(18) 



dfjr) 
dr 



and /'(r) 

At this point, we minimize the energy functional with respect to the two-body correlation 
function, /(r), to obtain the following Euler-Lagrange differential equation 



[^f{r)V{r) + 2Xf{r) 



l/7(r)^ 2- 
2V p 



d_ 

dr 



2f'{r) 



1 



1 /^x2nn 
2V p 



0, 



(19) 



where A is the Lagrange multiplier which imposed the normalization condition (\&|^&) = 1. 
By solving Eq. (fT9l) . using the numerical technique, the two-body correlation function, /(r) 
and therefore the energy of the system are obtained. This finally leads to the free energy 
function of the system 

F = E-TS, (20) 

where T and S are the temperature and entropy per particle of the systems (Fetter and 
Walecka, 1971). 

To calculate the free energy of liquid ^He, we use the Aziz interatomic potential (Aziz et 
al., 1979) in Eqs. ^ and ^ 



+ C-_ 
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10 



F{r) 



where 



F(r) 



^> D, 



(21) 



(22) 



and 



10.8K, A = 0.5448504 x 10^ 



a = 13.353384, 
Ce = 1.37732412, 
Cio = 0.178100, 



r„ = 2.9673A, 
Cs = 0.4253785, 
D = 1.241314- 



(23) 



The realistic Aziz Potential agrees with the He-He scattering experimental data which sat- 
isfies the following criteria: 



11 




FIG. 7: The free energy of liquid ^He as a function of number density at different temperatures. 

• It has a short-range repulsive part which described by exponential form 

• It has also a long range attractive tail includes the multipole interactions. 

A realistic potential between Helium atoms must have the criteria. Our results for the free 
energy calculations of the liquid ^He are given in Fig. ([7]) (Bordbar and Hashemi, 2002). 
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